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Lattice formulations of QCD with Wilson fermions and a chirally twisted quark mass matrix provide an at- 
tractive framework for non-perturbative numerical studies. Owing to reparameterization invariance, the limiting 
continuum theory is just QCD. No spurious quark zero modes, which are responsible for the problem with excep- 
tional configurations, can occur at finite values of the quark mass. Moreover, the details of the lattice formulation 
can be adjusted so as to simplify the renormalization and the 0(a) improvement of several quantities of phe- 
nomenological relevance. The first exploratory studies in the quenched approximation yield very encouraging 
results. 



1. Introduction 

In recent years a lot of progress has been 
achieved about lattice regularizations of gauge 
theories with fermions [ On one hand, local 
(but non-ultralocal) fermionic actions that enjoy 
a lattice form of chiral symmetry and (almost) 
ideal renormalization properties have been dis- 
covered and put at work [ for all of them 
the critical Dirac operator satisfies the celebrated 
Ginsparg- Wilson relation [ |[. On the other 
hand, more traditional formulations of lattice 
QCD (LQCD) based on Wilson and staggered 
quarks have been refined and widely used in re- 
alistic computations of hadronic observables and 
matrix elements. The simplest of these compu- 
tations are currently performed with dynamical 
quarks. 

In this contribution, I report on recent devel- 
opments about Wilson fermions, which remove 
practical obstructions in dealing with light quarks 
and simplify the renormalization and 0(a) im- 
provement of a number of observables. Improve- 
ments in the formulations of LQCD with stag- 
gered quarks are discussed elsewhere [ || . 

After shortly recalling the status of LQCD with 
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Wilson quarks and its problem with exceptional 
configurations (sect. 2), I introduce -for the case 
of Nf = 2 flavours- the formulation with chi- 
rally twisted mass (sect. 3). Then I discuss its ba- 
sic properties (sect. 4), the first non-perturbative 
studies (sect. 5) and analogous lattice formula- 
tions of QCD with Nf > 2 quark flavours, which 
can simplify the renormalization of some opera- 
tors of the effective weak Hamiltonian (sect. 6). 

2. Wilson fermions 

The well known lattice formulation introduced 
by Wilson [ || provides a gauge invariant reg- 
ularization for QCD with any number of quark 
flavours: the action is ultralocal and respects the 
global flavour symmetries of QCD, but all axial 
symmetries are broken by the Wilson term. This 
is no principle problem, as the flavour chiral in- 
variance SU(Nf) (g) SU (Nf) can be restored in the 
(quantum) continuum limit [ while the axial 
U(l) invariance of the classical continuum theory 
is broken by quantum fluctuations. 

The lack of chiral symmetry entails however 
some practically important consequences. First, 
complicated patterns of operator mixings arise, so 
that in many cases several operator subtractions^] 

2 The corresponding coefficients can be computed in per- 
turbation theory and beyond, e.g. by requiring the chiral 
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are needed in order to restore the chiral multiplet 
structure. Then, the leading deviations of renor- 
malized quantities from their continuum limit val- 
ues are linear in the lattice spacing a and typically 
non-negligible. Last but not least, the absence of 
a lower bound on the norm of the eigenvalues of 
the massive Dirac matrix may lead to unphysi- 
cal divergences in fermionic observables on non- 
trivial gauge backgrounds. 

2.1. O(a) improved Wilson fermions 

The problem with the leading cutoff effects has 
found a clean solution via Symanzik's improve- 
ment programme [ |, ||, , which allows to 
define and compute on-shell renormalizcd quanti- 
ties with leading cutoff effects of order a 2 , though 
at the price of tuning the coefficients of further 
countertermsQ. The LQCD action for the 0(a) 
improved theory with N{ quark flavours reads: 

S S [U; gl] + a 4 ]T fa) [PwM + M M (*) M 

X 

where S g [U; g$} is the pure gauge action with cou- 
pling M q = diag(mq, niq, mi, . . .) is the sub- 
tracted (M q = Mo — m c ) quark mass matrix and 

Dw[U] = {7V - §V*V + c sw f iaP} [U]+m c (2) 

represents the critical Wilson-Dirac operator. 
Omitting Lorentz indices, we denote by V (V*, 
V) the forward (backward, symmetrized) covari- 
ant lattice derivative and by the clover lattice 
discretization of F^. The coefficients c sw and 
am c are dimensionless functions^] of g^ and iVf. 

2.2. Exceptional configurations 

In quenched simulations of LQCD on relatively 
coarse lattices and for moderately light quark 
flavours, gauge configurations can be sampled 
where the Wilson-Dirac matrix has one or few 
eigenvalues with norm exceptionally close to zero, 

Ward identities to hold up to cutoff effects [ ^] . 
3 Implementing this programme may be non-trivial and 
CPU-time demanding, e.g. in the theory with non- 
degenerate and/or dynamical quarks, as well as for op- 
erators with complicated mixings [ h2[ . 
4 At the non-perturbative level, c sw and am c are uniq uely 
defined only up to O(a) and 0(a 3 ) corrections [ nil, re- 
spectively. 



i.e. much smaller than on the other configu- 
rations. The corresponding eigenvectors of the 
Wilson-Dirac matrix are referred to as "spurious" 
quark zero modes, because in a chiral invariant 
formulation of QCD the Dirac matrix can have 
zero modes only if some quark flavour is mass- 
less. 

On gauge configurations with spurious quark 
zero modes fermionic observables may undergo 
fluctuations that exceed in modulus the typical 
ones by orders of magnitudes: in this sense, these 
configurations appear to be "exceptional" . More- 
over, increasing the statistics does not reduce 
in general the standard deviation of the observ- 
ables, as further, even larger fluctuations may oc- 
cur. A reliable statistical analysis of the fermionic 
observable becomes hence impossible [ [ll], 
When employing the non-perturbatively 0(a) im- 
proved [ [ll) Wilson action, eq. (0), and the pla- 
quette gauge action, this problem is felt at values 
of m q of about one half the value of the strange 
quark's mass if a ~ 0.1 fm [ |l4). With the unim- 
proved Wilson action (c sw = 0) quenched compu- 
tations can be pushed to somewhat lower values 
of the quark mass, but cutoff effects are harder to 
control [ |l5| . 

In unquenched simulations, the importance 
sampling must give spurious quark zero modes 
a vanishingly small probability. State-of-the-art 
algorithms, however, account for the effects of 
sea quarks in a stochastic way: it is hence pos- 
sible that configurations with nearly zero-modes 
of the Dirac matrix are produced in the updating 
process (e.g. a HMC trajectory). This would 
then result in an exceptionally low acceptance 
rate(| which can possibly be avoided at the price 
of slowing down the algorithm. In the case of 
partially (un)quenched simulations, the problem 
with spurious valence quark zero modes is only al- 
leviated: exceptional fluctuations of hadron cor- 
relators have indeed been observed [ [l6) . 

A typical example of exceptional configurations 
is given in Fig. | There we plot versus the num- 
ber N of independent measurements the relative 

5 Exceptionally large values of the product 'driving force 
times MD integration step" have been observed [ Il7j| in 
simulations with the HMC algorithm. The eigenvalues of 
the Dirac matrix were not simultaneously evaluated. 
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Figure 1. Normalized standard deviation 

V~N(r(f) // vs. the number of measurements N in 
simulations with standard (empty symbols) and 
chirally twisted (filled symbols) Wilson quarks: / 
is a pion correlator at time separation t = 24a. 



standard deviation, multiplied by N 1 / 2 , of a pion 
channel correlator at fixed time separation: for 
the statistical analysis to be reliable, this quan- 
tity should approach a constant value as N ^> 1. 
Our example refers to a quenched simulation at 
13 = 6 (a ~ 0.1 fm) on a 24 3 • 48 lattice, with de- 
generate quark masses such that mpg/my — 0.47: 
employing the non-perturbatively 0(a) improved 
action ([[]) leads to untolerably large and rare fluc- 
tuations of the observable (empty symbols) . On 
the other hand, one can work at equivalent pa- 
rameters in a formulation based on Wilson quarks 
with chirally twisted mass, see sect. 3: of course, 
the bare quark mass parameters have to be ad- 
justed so as to match the value of the renormal- 
ized current quark mass (to^ in the notation of 
sect. 4). In this case the occurrence of spurious 
quark zero modes is avoided: as shown in Fig. [I] 
(filled symbols), there are no exceptional fluctua- 
tions in the MC history of the pion channel corre- 
lator, so that a reliable statistical result^] can be 
obtained. 

3. Chirally twisted "mass" 

The presence of chirally twisted "mass" terms 
is a general feature of lattice gauge theories 
with Wilson fermions, as soon as the "physical" 
fermion mass term is not aligned in chiral space 
to the Wilson term. Of course, the latter gener- 
ates itself a "mass" term, whose divergent part 
requires appropriate subtraction. 



3.1. Twisted mass and parity 

To the best of my knowledge, Osterwalder and 
Seiler first proposed [ [l!| a lattice formulation 
where the fermionic action is of the form 

S^=^'(x) [{D%+fhW] (x), 

D% = 1 V + e-^[mc - § V*V + c sw fiaF] , (3) 

namely with <f> = n/2 and csw = 0. The angle <j> 
can be taken as a measurement of the misalign- 
ment between the fermion mass term (oc in) and 
the Wilson term, subtracted by a critical mass 
and a clover-like termQ It was also observed that 
lattice formulations with fermionic action (|3|) in 
general do not respect parity [ ^| . Indeed, as long 
as (f> 7^ (mod tt), quantum fluctuations generate 
in the effective action a term oc tr (FF), which 
survives to the continuum limit (if not canceled 
by a suitable counterterm) . 

Parity breaking is an unwanted feature if one is 
interested in QCD. One can show that, up to cut- 
off effects, parity breaking is avoided by consid- 
ering Nf > 1 quark flavours and giving the array 
4> = . . . , 4>( Ni ^] a non-trivial flavour depen- 

dence such that X)/=i 4^ — 0- ^ n the case of 
N{ — 2, which is relevant e.g. for u and d quarks, 
one is then led to the fermionic action: 



Si = a 4 l'( x ) ( D tmw + m q) <l'( x ) ( 4 ) 



where q' is a quark doublet field, the quark mass 
term (oc m') takes the usual form, but the critical 
Wilson-Dirac matrix is chirally twisted^]: 

DtmW = 7V + e-^ 3 [m c - f V* V + c sw f iaF] (5) 
A (non-anomalous) change of quark field basis, 



q = q> e -*^ 3 / 2 , ( 6 ) 



3 The full statistics was of N = 260 measurements 



brings the action (Q) into the form 

Si = a 4 Y q(x) (-t>w + m q + *Mq75T 3 ) q(x) (7) 



where fi q = m' q sin u> is the twisted mass parame- 
ter and m q = nriQ—m c = m' q cosw. The bare mass 

7 In general, am c and cgw may differ from am c and c sw . 
8 The matrix t 3 in fJr mW acts in flavour space. 
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parameters are /i q and m , as the m c -dependence 
cancels between m q and D^. 

The formulation with fermionic action (Q) 
or ^ is referred to as lattice twisted mass QCD 
(LtmQCD). At the classical level, it obviously 
represents a regularization of QCD with Nf = 2 
(mass-degenerate) flavours. This property holds 
true at the quantum level [ see sect. 4. 

3.2. Protection against exceptionals 

The lattice theory with fermionic action (Q) has 
been studied (for c sw = 0) to establish whether, 
at finite g^ and for values of aing in a certain 
range, parity and isospin can be spontaneously 
broken^, as conjectured by Aoki [ ^2j. It was 
noted [ |2^] that the lattice Dirac matrix corre- 
sponding to eq. (0) can not be singular on any 
gauge background as long as /i q ^ 0, since 



< Bct\D 



w 



^751 



det[(D^ + m,)*^ + m q ) + fx 



(8) 



Here Det[. . .] (det[. . .]) denotes the fermionic de- 
terminant in the two-flavour (one-flavour) space. 

The one-flavour Dirac matrix corresponding to 
the action (|J) is also not singular on any gauge 
configurations, as long as rh sin <\> ^ 0. This is be- 
cause |det[Z?^ + 77T.] | 2 can be written in the same 
form as the r.h.s. of eq. (||) with m q =>• rhcoscj) 
and fj, q =>■ rh sin <j). Based on this property, 
the authors of Ref. [ Q proposed to employ 
the fermionic lattice action (^) with rh c = m c , 
csw = c sw , ir/2 > \(p\ > and m > to avoid 
the occurrence of spurious quark zero modes in 
massive LQCD with Wilson quarks. They pre- 



sented numerical evidence, see Fig. 3 of Ref [ 24 



that exceptional configurations are avoided in the 
quenched approximation and discussed on a semi- 
classical level the relation between their lattice 
formulation and QCD in the continuum limit. It 
should be noted that in the quenched model the 
parity breaking inherent to the formulation with 
action (||) reduces to a mere 0(a) effect on renor- 
malized quantities, but this is no longer true with 



9 In this context /x q played the role of an external zero- 
momentum source, which can orient the vacuum in a cer- 
tain direction (in flavour chiral space), and the main inter- 
est was on the properties of the vacuum after taking the 
thermodynamic and zero-source limits. 



(partially) unquenched quarks. 

4. Basic properties of lattice tmQCD 

I discuss here the properties of LtmQCD as an 
ultraviolet regularization of QCD with Nf = 2 
mass-degenerate quark flavours^]. 

4.1. Symmetries and renormalizability 

The LtmQCD action, see eq. (||) or (]?]), is in- 
variant under lattice gauge transformations and 
translations, axis permutations, charge conjuga- 
tion and the global symmetry U(l) correspond- 
ing to conservation of the fermionic number. At 
/Lt q ^ <^> u> — &rct&n(ii q /m q ) ^ 0, isospin sym- 
metry is reduced to a U(l)-invariance with gen- 
erator r 3 /2 and axis reflections, such as parity, 
are no longer a symmetry. It is important to note 
the residual invariance Pp, 

U (x) -» U {x P ) , U k {x) -» U^ 1 {x P - ak) , 
q(x) -> T 1 7 q(a;p) , q{x) -> q{x P )^ T 1 , (9) 

where xp = (xo, — x), as this symmetry rules out 
a term cx tr(FF) in the quantum effective action. 

A standard analysis [ |2lJ based on lattice sym- 
metries and power counting shows that the model 
is renormalizable. The relations between bare 
and renormalized action parameters take in gen- 
eral the form (I recall m q = mo — rn c ): 



Mr 

TOR 



= ^g(ffo, a TO q , aAi q ; aii) gl , 
= Z^(gl,am q , aii q ; a/i) Mq , 
= Z m {gl , am q , a^ q ;aii)m. 



q ' 



(10) 



but the renormalization factors Z can be chosen 
to be independent of am q and api q (mass inde- 
pendent schemes). The additive renormalization 
of Too is independent of the quark mass param- 
eters f @, |25|, |2l: m 



|25| , |26j| : m c = m c (gl), up to intrinsic 
O(a) corrections in the case of non-perturbative 
determinations of m c . The ratio Z^/Z^ does not 
depend on the subtraction scale [i, consistently 
with recovery of flavour chiral symmetry in the 
continuum limit. Noting the bare lattice identity 



d*X = -2Mqe 36c P 



(11) 



10 A satisfactory regularization of Nf = 2 QCD with bare 
quark masses m' q ± Sm' q can be obtained e.g. by adding 
to the action density in eq. (Q) a term q' (x)t 1 Sm' q q' (x) . 
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where V b denotes the one-point split vector cur- 



1,2) 



rent, it is easy to argue [ 21 that Z^Z R = 1 



4.2. Continuum limit and cutoff effects 

Although the physical interpretation of the 
fermionic correlation functions is most transpar- 
ent in the quark basis corresponding to eq. (^), 
the renormalization of gauge-invariant correla- 
tion functions, including those with insertions of 
local operators, looks simpler in the quark basis 
corresponding to eq. (fjj). In this basis, indeed, 
the critical Wilson-Dirac matrix is given by 
and the mixing properties of the operators in the 
massless theory (fi q = 0, m — m c ) are usually 
well known: the construction of lattice fields that 
are multiplicatively renormalizable and respect, 
up to cutoff effects, the chiral multiplet structure 
is hence straightforward. Concerning the overall, 
possibly scale-dependent renormalization factors 
of the various operator multiplets, we assume for 
simplicity [ ^| that they, as well as Z & , Z m and 
Zfj,, are chosen to be independent of 10. 

In the quark basis of choice, the Ward identities 
of flavour chiral symmetry read (for b — 1,2,3): 



2m R (P R ) b + z^ R 6 3b (S R )° 



d^(V R )l ~ -2 MR e 3bc (P R ) c , (12) 

where the symbol ~ means equality up to 0(a) 
corrections and renormalized operators are de- 
fined as usual with Wilson fermionspl, e.g. 



( {V R , A r } t = q{Z vlfl , Z A%ll5 }± 



t q. 



(13) 



The form of eq. ( ^2[ ) reminds us that appro- 
priate linear combinations of renormalized opera- 
tors, reflecting the formal change of variables (||) , 
should be taken, in order to obtain in the con- 
tinuum limit operators with well defined parity 
and isospin properties. Exploiting eq. ( |l2| ) one 
can indeed extract m R and fi R and introduce the 
renormalized counterparts of m' q and lu: 



m R + Mr ' a = arctan 



\m R J 



(14) 



I remark that tan(a) = {Z^IZ m ) tan(w) and 



rrv x 



z m >(g 2 



Ttq. Then, if one considers (with 

11 For 6 = 3 severe power-like divergences must be sub- 
tracted to define which is no problem in perturba- 
tion theory but is delicate at the non-perturbative level. 



(A' R t 



(Vk)1, (A^l = (A R )l 



l B.JfJ, 



(pkr 



the renormalized operators 

cos(a)(V R )l 
cos(a)(A R )l 
cos(a)(P R f 



, sin(a)(A R )^ 
e 3ic sm(a)(V R )l 
§sin(a)(S R ) 



(15) 



the Ward identities of flavour chiral symmetry 
take the usual form [ Ell : 



d,(A R 



M - R , M -r 2m' R (P R ) b , dM% ~ . (16) 

This result reflects the existence of a linear 
mapping among renormalized correlation func- 
tions computed with action parameters that cor- 
respond to different values of a and identical val- 

The mapping between 



ues of g R and m! x 



R 
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correlators of gauge invariant fields with unequal 
space-time arguments at a = a and a = 
readspl: 



(tf (*)•••>, 



.0 ! 



>(x) = R%>(a)4>£\x). (17) 
(r) 

Here (jr k ^ denotes the /c-th field component of the 
renormalized chiral multiplet r and R( r \u>) is the 
(formal) multiplet transformation matrix under 
quark transformations of the type The dots 
in ( |l7| ) stand for further local fields, with those 
on the l.h.s. being related to those on the r.h.s. by 
the appropriate product of chiral transformation 
matrices R. The relations ( [l7| ) are regularization- 
independent properties and imply the recovery 
of parity, isospin and chirality in the continuum 
limit of LtmQCD [ |2^|. Moreover, it is useful to 
note that they look just as one would guess from 
simple formal arguments. 

The full information on QCD with two mass- 
degenerate quark flavours can hence be obtained 
by computing correlation functions with the ac- 
tion (jij) or (0). It is understood that the contin- 
uum limit must be approached at fixed values of 
the renormalized parameters, including a. In this 
way, the cutoff effects inherent to any determina- 
tion of a at finite lattice spacing have no impact 
on the (extrapolated) continuum limit results^ 

12 The definitions (|l5|) are special cases of those in eq. ([it]), 
13 Phenomena like a possible Aoki phase [ 22l Gq], are im- 
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4.2.1. Simplified operator mixings 

The regularization of QCD with fermionic ac- 
tion (^) differs from Wilson's (improved) one as 
long as the quarks are massive. In the massless 
theory, fi q — and tuq = m c , any difference dis- 
appears, implying that the pattern of leading op- 
erator mixings is globally the same as in the orig- 
inal formulation by Wilson. 

However, given a certain fermionic field, its 
physical interpretation depends on how it trans- 
forms under the parity and isospin operations, 
whose form is in turn dictated by the parame- 
terization of the quark mass term. On the other 
hand, the mixing properties of the fermionic field 
depend crucially on its chiral orientation with re- 
spect to the Wilson term. Hence, for a specific 
fermionic field with a given physical interpreta- 
tion (e.g. the physical axial current), the mixing 
properties may change, even at the leading level, 
with us (or a). For instance, eq. ( |l5| ) shows that 
the component (A'^) 1 ^ of the physical non-singlet 
axial current is given by (Ar)* if a = lu = and 
by {Vr)1 if a = lu = tt/2. The currents (A R )l 
and (Vr)„, which are defined in eq. ([l3|), involve 
different renormalization factors Za and Zy, with 
Zy = 1 if the one-point split definition of the bare 
current V? is adopted. 

One can then adjust the ultraviolet regular- 
ization so as to simplify as far as possible the 
renormalization properties of certain -but not 
all- physical observables of interest. E.g. work- 
ing with lu = 7r/2 entails remarkable simplifica- 
tions [ in the renormalization of the decay 
constant of the charged pseudoscalar mesons and 
the chiral condensate, as (^4r)^ = (Vr)^ and 
(S R )° = 2z(F R ) 3 . 

4.2.2. Hamiltonian formalism at fixed a 

Unimproved (c sw = 0) lattice tmQCD admits a 
positive and selfadjoint transfer matrix [|2(J, with 
the constraint |8+2amo| > 6 and no constraint on 
a/Xq. Lattice correlation functions can hence be 
represented as usual in terms of operator matrix 
elements with time-dependent coefficients. 

From the symmetry properties of LtmQCD, see 

material for any ^ 0, if the continuum limit of renor- 
malized quantities is taken at fixed renormalized parame- 
ters. 



sect. 4.1, it follows that the set of lattice quan- 
tum numbers for lu ^ is reduced as compared 
to that of Wilson's formulation (lu — 0): parity 
and isospin are replaced by the quantum num- 
bers pf and q\ corresponding to the unphysical 
parity transformation Pp, eq. @, and the un- 
broken isospin generator. This implies, for in- 
stance, that in LtmQCD the physical vacuum and 
neutral pion states are labelled by the same set 
of lattice quantum numbers and can be interpo- 
lated by the same lattice field. Moreover, lattice 
operator mixings are constrained by pp and q\ 
rather than parity and isospin quantum numbers. 
For instance, the operator (^4r)^, which yields, 
as a — > 0, the first component of the physical 
isotriplet axial current, mixes at order cs/zr with 
the second component of the physical isotriplet 
vector current. 

In the quantum mechanical analysis of tmQCD 
correlators (lu ^ 0) at fixed a, one must hence 
include [ |T[ [l^] the contributions of matrix ele- 
ments that would vanish if parity and isospin were 
exact symmetries)^]. The 0(a)-improvement is 
expected to reduce the size of these contributions. 

4.3. O(a) improvement: the case a = tt/2 

The O(a) improvement of LtmQCD [ |f| |^ 
is most conveniently discussed in the unphysical 
quark basis corresponding to eq. (0). Without 
loss of generality, one can assume that an infrared 
cutoff is in place (thanks to some specific choice 
of external momenta or boundary conditions), so 
that all correlation functions admit a Taylor ex- 
pansion in powers of /x q and m q around the mass- 
less theory. Once the latter has been fully renor- 
malized and O(a) improved in a mass indepen- 
dent scheme, the renormalized action parameters 
of the massive theory are consistently [ [n| de- 
fined by eq. (|l0|), with g^, m q and /i q replaced 
by: 

ffo = .9o(! + h am <\) ) 
m q = m q + b m am q + b m a/i q , 

= Mq(l + b fl am q ) (18) 

14 The renormalized counterparts of these matrix elements 
actually vanish like a^tR, or faster, in the continuum limit, 
because of the relations (II 71). 
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and Z s , Z m and Z M depending only on g 2 . The 
absence of a term oc a^i q in the expression for 
5q reflects the fact that the partition function 
of LtmQCD is even in pb q , see eq. for the 
fermionic determinant: hence an action countert- 
erm of the form afj, q tT(FF) can not be generated. 

Improved operators can be obtained by appro- 
priate subtraction of the operator mixings that 
come with powers of am q or a/i q . In absence of 
power divergent mixings, the construction of mul- 
tiplicatively renormalizable improved operators is 
rather simple, e.g. 



(Pit 



an q b A e 3b V^ 



= P 



b = 1,2. 



(19) 



The corresponding renormalized and improved 
operators are obtained after rescaling by suitable 
factors Za(1 + bp^arn^ and .Zp(l + 6pam q ). It is 
useful to remark that the set of 0(a) counterterms 
that have been introduced is slightly redundant [ 
p6[ : one of them can be chosen arbitrarily, e.g. 
by setting b m = —1/2. 

Adopting a formulation with |m q | ~ |m q | < 
0(a) and fi q ^ 0, i.e. 



7r/2 + 0(a). 



(20) 



we find that <?j| = Z g g^ and = Z^fi q up to 
0(a 2 ) corrections, see eq. (|l8|), and all 6-type im- 
provement coefficients are not needed. One needs 
instead to know the 6-type improvement coeffi- 
cients, like oa in eq. ([lj]), which are associated to 
operator mixings that violate parity or isospin. 
This observation implies that one can determine 
the 6-type coefficients by requiring the improved 
correlation functions to satisfy parity and isospin 
invariancesp[ Implementing the condition (20) 
requires to know only c sw and m c (unimproved). 
Moreover, for a as in eq. (|20|), eq. ([ll] ) is inter- 
preted as the physical PC AC relation, see eq. (p^|) : 
O(a) improved estimates of the mass, mpg, and 
the decay constant, Fps, of the charged pseu- 



15 In principle, it is also possible to get rid of the mix- 
ings of order ap,^ without improving the operators: this 
is the case if in the quantum mechanical analysis of the 
correlators one can disentangle, e.g. by looking at the val- 
ues of the effective masses, the contributions arising from 
intermediate states with different parity and isospin. 



doscalar meson can hence be obtained [ ^7J from 
the correlator 

/i q G P (xo) =a 3 ^ / i q (P 1 ( a; )P 1 (0)), (21) 

X 

if c sw and m c are known. If one wishes to work 
with | a | = 7r/2 + 0(a 2 ), an improved estimate of 
m c (which may require to know ca) is also needed. 

5. The first non-perturbative studies 

The above theoretical understanding of 
LtmQCD has been checked at the one-loop 
level [ |2f|, as far as the action and the oper- 
ators of eq. ( |l5| ) are concerned. At the non- 
perturbative level, test studies have been carried 
out in the quenched approximation]^ using the 
non-perturbatively improved [ |ll| Wilson action. 
The authors of Ref. [ ^ studied LtmQCD in 
a box of volume L 3 x T with Schrodinger func- 
tional boundary conditions [ |2fi|] , L = 0.75 fm, 
T = 2i, Lm' R = 0.154, a = 1.44 and four lattice 
resolutions: a € [0.093,0.046] fm. The scaling 
behaviour of some renormalized and improved 
quantities (which in large volume yield the mass 
and the decay constant of pseudoscalar and vec- 
tor mesons) was found to be consistent with O(a) 
improvement, with the residual cutoff effects at 
a = 0.093 fm ranging from 0.5% to 9%. 

The same setup and observables have then been 
employed for a study in realistically large vol- 
umes [[ED: L = 1.5 to 2.2 fm and T = (2 3)i, 
so to ensure mpsL > 4.5. This study was re- 
stricted to two lattice resolutions, a — 0.093 
and 0.068 fm, and, for each of them, three sets 
of quark mass parameters, which correspond to 
|a| — n/2 + 0(a) and pseudoscalar meson masses 
in the range 1.85 > (mps/wK±) 2 ^ 0.85. 

Another test study of quenched LtmQCD with 
non-perturbatively improved action has been per- 
formed at fixed spatial volume, L = 1.5 fm, with 
T = 2L, lattice spacing a = 0.093 fm, peri- 
odic boundary conditions and several values of 
the quark mass [ ^TJ : the lowest value of mpg is 
about 320 MeV, though with m PS L ~ 2.4. In 
none of these studies exceptional configurations 



ie For all quenched studies I conventionally set the "phys- 
ical" scale by employing Sommer's scale: rg = 0.5 fm. 
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• Ref. [18] 
> Ref. [30] 



Figure 2. r rn^, s /m' K vs. r m' R at j3 = 6. The 
dashed vertical line corresponds to mpg = mx± . 

were found, while values of mps well below TOk± 
were reached with a moderate computational ef- 
fort: results for Ttps/ m R at a = 0.093 fm from 
Refs. [ H are shown in Fig. |. 



0.38 
-0.01 







• Ref. [18], L=1.5fm 






Ref. [18], L=2.2fm 






:: Ref. [14], L=1.5fm ~ 












-_f-- 






I 



0.02 

(a / O 2 



Figure 3. Fps vs. a 2 at tops — 1.2mK±- The con- 
tinuum extrapolation by Ref. [ Q is also shown. 

The scaling behaviour of Fps in large volume is 
presented in Fig. || for the O(a) improved Wilson 
formulations with \a\ ~ tt/2 [ |l§ and a = [[uj. 
In the latter case, where four lattice resolutions 
were considered to allow for continuum extrapo- 
lation, the simulation data have been reanalysed 
to precisely match the renormalization conditions 
adopted for LtmQCD. Following closely Ref. [|l4|, 
an analogous comparison has been carried out, 
see Fig. 0, for the combination of renormalization 
group invariant quark masses M + M s , where M 
is the average mass of the u and d quarks. The 
results for Fps r and (M + M s ) r that are ob- 
tained from the two lattice formulations should 
agree in the continuum limit: this seems to be 
the case within the statistical errors shown in the 
figures^. Moreover, in agreement with the in- 
dications of the scaling test at L ~ 0.75 fm, the 
tmQCD estimates of these quantities show pretty 
small cutoff effects. 



17 For LtmQCD m^/rripg at the point mpg = m K ± could 
be obtained by simple interpolation of simulation data. 



0.38 
0.38 - 
0.34 
0.32 



• tmQCD, /S=6.2- 
. tmQCD, (3 = 6.0 



i stQCD 



3 stQCD-CL 



[NPB571(2000)237] 



0.02 
(a/r ) 2 



Figure 4. (M + M s ) r vs. a 2 ' jr 2 for a = 



(stQCD: [hj) and \a\ ~ tt/2 (tmQCD: [hj). 



Standard linear solvers (CGNE, BiCGStab) 
were employed to invert the (preconditioned) 
LtmQCD Dirac matrix. CGNE worked fine, 
while BiCGStab often got stuck for \a\ ~ 
| a; | ~ 7r/2. The implementation of the 
twisted mass term /z q in unquenched algorithms 
(HMC, PHMC, GHMC, multiboson) is straight- 
forward^]. 

6. B K and the AI = 1/2 rule 

The two-flavour formulation -LtmQCD- that 
I discussed above can be extended in a vari- 
ety of ways to describe QCD with several non- 
degenerate flavours of Wilson quarks. The de- 
tails of the (Wilson-like) fermionic regularization 
should be fixed, case by case, so to render as sim- 
ple as possible the mixing pattern of the operators 
that are relevant for the physical problem of in- 
terest, while avoiding the occurrence of spurious 
quark zero modes. The general reasons underly- 
ing this possibility are discussed in sect. 4.2.1. 

For the computation of Bk (K°-K° mixing) 
it is convenient to adopt the lattice formulation 
with fermionic action density: 

C F = q(D^ + im' qll5 T 3 )q + s(D^ + m' qs )s (22) 

which corresponds to take a = tt/2 for the two de- 
generate light flavours and a — for the strange 
quark. In this quark basis the operator for (phys- 
ical) parity-even AS — 2 transitions reads [ [2l| | 

O'jyTA = -Hsi^d){- Sllll5 d) = -2W sd A sd {22.) 



18 Unquenched simulations with the action (|(j) have been 
performed to study issues related to the Aoki phase [ . 
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and is multiplicatively renormalizable. A compu- 
tation of i?K based on this approach is currently 
in progress [ ^] . 

For the computation of CP conserving AS = 1 
matrix elements of the weak effective Hamilto- 
nian with active charm flavour, it is convenient 
to adopt the following lattice formulation [ [33) : 

C F = $ (D'w + m q + i£ q j 5 T 3 ) ip (24) 

with vp T — (u,d; s,c) and quark mass matrices 

m q = diag(m ql ,m q i;ra qs ,TO qc ) 

/iq = diag(/i q l, -Mq^Mqs^qc) • (25) 

Thinking of -0 as a pair of doublets and setting 

—3^- = tan uji , qs = — —3-2- = tancj/j , (26) 

m ql m qs m qc 

one can show [ [33) -by arguments analogous to 
those of sect. 4.2- that the operator 

°VA+AV = {{VsdAuu + A sd V uu ±...)-{u^c)] 

must be interpreted as the parity-even operator 
i 0' V y +AA (relevant for K — ► it transitions), if 
u>i = u>h = 7r/2, and as the parity-odd operator 
O'va+av (relevant for K — > irir transitions), if 
loi = —uJh = 7r/2. Symmetries and power count- 
ing imply that Oy A+AV does not mix with di- 
mension six operators, while the leading mixing 
with dimension three operators has a coefficient 
oc a (fXq C — jUqi)(/iq S — Mql)- Quadratic diver- 
gences in K — > 7r matrix elements and spurious 
quark zero modes in general are hence avoided. 

7. Conclusions 

Wilson fermions with chirally twisted mass can 
provide a variety of alternative regularizations for 
lattice QCD. Thanks to an up to now unexploited 
freedom in formulating lattice QCD with Wilson 
quarks, they allow to avoid some of the most se- 
rious problems due to lack of lattice chiral sym- 
metry, while preserving ultralocality of the action 
and hence a moderate computational cost. 
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